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Abstract. A fully self-contained model of homochirality is presented that contains the effects 
of both polymerization and dissociation. The dissociation fragments are assumed to replenish 
the substrate from which new monomers can grow and undergo new polymerization. The 
mean length of isotactic polymers is found to grow slowly with the normalized total number of 
corresponding building blocks. Alternatively, if one assumes that the dissociation fragments 
themselves can polymerize further, then this corresponds to a strong source of short poly- 
mers, and an unrealistically short average length of only 3. By contrast, without dissociation, 
isotactic polymers becomes infinitely long. 
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1. Introduction 

Central to the question of the origin of life is the polymerization of the first 
complex molecules that can have catalytic properties and that would even- 
tually carry genetic information. It is widely accepted that our current life 
form involving DNA carrying the genetic code and RNA producing the pro- 
teins that, in turn, catalyze the production of nucleotides, must have been 
preceded by a simpler life form called the RNA world (Woese, 1967; Crick, 
1968; Orgel, 1968; see also Wattis & Coveney 1999). Here, the RNA has 
multiple functionality, it carries genetic code and it is also able to catalyze 
the production of new nucleotides. 

The RNA of all terrestrial life forms involves a backbone of dextrorotatory 
(right-handed) ribose sugars. Theoretically, life could have been equally well 
based on levorotatory (left-handed) sugars. Unless this selection was some- 
how externally imposed, e.g. via circularly polarized light (Bailey, 2001), 
magnetic fields (Thiemann 1984), or via effects involving the parity-breaking 
electroweak force (e.g., Hegstrom, 1984), this must have been the result of 
some bifurcation process. Indeed, the homochirality of left-handed amino 
acids and of right-handed sugars in living cells can be explained as the re- 
sult of two combined effects, auto-catalytic production of similar nucleotides 
during their first polymerization events and a competition between left- and 
right-handed nucleotides. The general idea goes back to early work of Frank 
(1953), and has been developed further by Kondepudi and Nelson (1984), 
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Goldanskii and Kuzmin (1989), Avetisov and Goldanskii (1993) and more 
recently by Saito and Hyuga (2004a). Of particular interest here is the re- 
cently proposed detailed polymerization model of Sandars (2003); see also 
Brandenburg et al. (2005, hereafter referred to as BAHN) and Wattis and 
Coveney (2005). The main point of Sandars' model is the assumption that 
the polymerization of monomers of opposite handedness terminates further 
growth on the corresponding end of the polymer. This is referred to as enan- 
tiomeric cross-inhibition. Such inhibition makes it generally quite hard for 
any polymer to grow successfully. However, once a polymer has become 
successful in reaching an appreciable length, it will have catalytic properties 
promoting the production of monomers of the same chirality as that of the 
catalyzing polymer. 

All the polymerization models presented so far ignore the possibility of 
polymers breaking at an arbitrary location. Without this process polymers 
would, in the homochiral case, grow to infinite length which is clearly un- 
realistic. We begin by discussing a model for the dissociation of isotactic 
polymers, where all the building blocks have the same chirality. Next, we 
consider the dissociation of polymers whose one end has already been spoiled 
with a monomer of the opposite chirality. We then incorporate the dissocia- 
tion model into the full polymerization model of Sandars (2003) and discuss 
an important modification that is necessary to prevent the average polymer 
length from being too short. 



2. Outline of the model 

The model that we are proposing has arisen through the realization that the 
obvious generalization of the polymerization model of Sandars (2003), to 
include dissociation, leads to two important difficulties. It was our desire to 
resolve these problems in a way that seemed most natural to us, and that 
involves the least amount of assumptions and new parameters. What we came 
up with is a closed model that is fully self-contained. As in the original 
model of Sandars, new monomers of either chirality are being produced from 
an achiral substrate. However, unlike the original model, no external source 
of the substrate is required. Instead, the substrate can be replenished by the 
"waste" generated by fragmented polymers. 

Before we can discuss the dissociation model, let us explain in a few 
words the polymerization model of Sandars. Here, polymers can grow by the 
addition of monomers that can have either the same or the opposite chirality, 
and the corresponding reaction coefficients are ks and kj, respectively. The 
subscript S indicates that the chirality of both reaction partners is the same. 
The addition of a monomer of opposite chirality leads to the inhibition of 
further growth at that end of the polymer, which is indicated by the subscript 
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/. The process of such an inhibition, also referred to as "enantiomeric cross- 
inhibition", is the single most important aspect of the model without which 
there would be no bifurcation from a racemic (i.e. equally many right and left 
handed building blocks) to a homochiral state. 

The fragmentation involves a new parameter: the decay rate 75, at which 
a polymer can break up anywhere in the chain. Again, the subscript S refers 
to the situation where the partners involved in the bond have the same chi- 
rality. If the chirality is different, we call the decay rate 77-, in analogy to the 
corresponding reaction coefficient kj in the original polymerization model of 
Sandars (2003). 

The perhaps most obvious assumption for dissociation would be to let 
the fragments continue to polymerize with new monomers. This leads to two 
undesired features of the model. In the previous case with only polymeriza- 
tion the homochiral equilibrium had the property that polymers of different 
lengths are all equally abundant. This goes on all the way to infinity. If we 
now allow these polymers to break, there is potentially a catastrophe in that 
arbitrarily many short polymers can form. This is also supported by the nu- 
merical simulations discussed below. Furthermore, the numerical solutions 
show that, even in the best possible case, the average polymer length never 
exceeds 3, which is clearly unrealistically short. We propose two alternative 
ways to allow for the formation of longer chains. One possibility is to in- 
clude an additional degradation of polymers leading to a loss term in the 
polymerization equations and a corresponding source term for the substrate. 
Another possibility is to recycle the dissociation products into the substrate 
without invoking an additional degradation of polymers. In both cases the 
total number of building blocks in the system is constant, so the substrate 
plays now an integral part of the model. As a mechanical analogue, we can 
think of the mass of the substrate as being similar to potential energy, and 
the mass of all polymers as being similar to kinetic energy, such that the total 
number (corresponding to the total energy) is conserved. Thus, not only goes 
the production of new left and right handed building blocks at the expense of 
the substrate, but now the substrate is being replenished by the dissociation 
fragments such that the total number of building blocks (regardless of their 
chirality) remains constant. 

In the following we develop the model step by step. We first need to show 
that the mean polymer length is never more than 3 if the fragments are reused 
for further polymerization. 



3. Developing the dissociation procedure 

Following the basic idea behind Sandars' polymerization model we assume 
the presence of left and right handed polymers of length n, denoted by L n 
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and R n , respectively. We also assume the presence of polymers whose one 
end has been spoiled by a reaction with a monomer of opposite chirality. The 
resulting polymers of this form are denoted by L n R\ and R n L\. 

3.1. ISOTACTIC DISSOCIATION 

We begin with the description of a dissociation model by discussing isotactic 
polymers of length n, which are assumed to break (dissociation) at a mean 
rate 75 (assumed independent of n), at position m, according to 

L n L m + L n _ m . (1) 

Here, L refers to left4ianded building blocks, but a corresponding equation 
is also valid for right4ianded polymers, denoted by R. In the present case, 
the fragments L m and L ra _ m will be reused for further polymerization. As 
an example, L4 can break up into two L2, or into one L\ and one L3, but 
for the latter there are two possibilities to do this. Thus, for n = 4 there are 
altogether n — 1 = 3 different ways of destroying L4. This then leads to an 
evolution equation for the concentration of polymers, [L n ], 

d N 

-[L„] = ... + 2 7s E [L m }-(n-lhs[Ln], (2) 

m=n+l 

where n > 2, and the last term represents the decrease of the concentration 
[L n ] due to the n — 1 different ways of breaking up the polymer. The first term 
represents the corresponding gain from breaking up polymers with m = n — 1 
or more building blocks. The evolution equation for [L{\ has only a gain term 
from breaking up polymers of length n > 2, so 

d N 

-[L 1 ] = ... + 2 75 ^[L n ]. (3) 

n=2 

The absence of any negative terms (sinks) on the right hand side implies that, 
if there is only dissociation, [L\] can only grow. The dots in Eqs. (2) and (3) 
denote the possible presence of extra terms (discussed in the next subsection) 
that would be needed to model the primary polymerization process. 

The same set of equations (2) and (3) applies also to R n . The mean rate of 
dissociation is again 75, so the model is completely symmetric with respect 
to exchanging L +± R. Using the identity 

N N N 

E [£n] = ££(n-l)n[L n ], (4) 

n=l m=n+l n=l 

one can easily see that these reaction equations (2) and (3), in the absence 
of extra terms, conserve the total number of left and right handed building 
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Figure 1. Evolution of [L n ], using as initial condition [Lioo] = 1 and [L n ] = for n ^ 0. 
blocks, i.e. 

N N 

El = n [L n ] = const, Er = ^ n[R n ] = const. (5) 

n=l n=l 

As an illustration we show in Figure 1 a numerical integration of the evolution 
of [L n ], using as initial condition [Lioo] = 1 an d [L n ] = for n ^ 100. Thus, 
we have El = 100 initially, and this value is preserved by the model for all 
times. 

As can be seen from Figure 1 , both monomers and short polymers are im- 
mediately being produced. For n > 2 the concentration reaches a maximum 
at a time that is of the order 7^ \ and decays then exponentially to zero. 

3.2. Semi- spoiled polymers 

For polymers whose one end has been spoiled by a monomer of opposite 
chirality, we have two types of reactions: those where the spoiling enantiomer 
breaks off (rate 7/) and those where the polymer breaks up somewhere else 
in the isotactic part (rate 75). Thus, we assume 

L n Ri L n + Ri, (6) 

and 

L n Ri — ► L m + L„_ m i?i (7) 

for 1 < m < n — 1. Ignoring a particular complication that will be discussed 
in a moment, our preliminary set of equations for these additional reactions 
is then given by 

d N 

^[LnRi] = - + 75 [L m Ri}-{li + (n-lhs}[L n Ri}, (8) 

m=n+l 
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d N 

— [L n ] = ... + 7S J] ^m^+T/^i], (9) 

m=n+l 

d 



dt 



[L 1 ] = ... + 7s ^M 1 ]+ 7; [L 1 fi 1 ], (10) 



n=2 

N 



^[R 1 ] = ... + 1I J2[L n R 1 ]. (11) 

n=l 

These equations ignore the dissociation of isotactic polymers discussed in the 
previous section, but they can simply be added to the present set of equations. 
Again, the system of equations has to be completely symmetric with respect 
to exchanging L +± R. However, the reaction (7) for m = n — 1 produces 
L\R\ at a rate that is proportional to [L n R\]. In general, since [L n Ri] / 
[R n Li], this would lead to [Lii?i] / [R\Li], which is not permitted. We 
therefore have to discard the reaction (7) for m = n — 1, i.e. we have to 
discard the reactions 

L n R 1 L n _i + L1R1 (discarded), (12) 

and likewise for the dissociation of R n L\. Since we have therefore one reac- 
tion less, this means that in Eq. (8), which now applies only for n > 2, the 
n — 1 factor changes effectively into an — 2 factor. Furthermore, in Eqs. (9) 
and (10) the sums start only with m = n + 2 and n = 3, respectively. 

When writing down the full set of equations we have to treat the evolution 
of [L1-R1] separately, so 

j t [L 1 R 1 ] = ...- 1I [L n R l l (13) 
while for n > 2 we have a pair of equations 

^[L n Ri] = - + wi L ^ - { 7/ + [n - 2) 7S } [LM, (14) 



^[RnLi] = ... + - { 7J + (n - 2) 7S } [Rn^]. (15) 



Here, 



w 



N N 

^= 7S E \ L ™Ri\, ^n RL) =7s E i R M (16) 

m=n+l m=n+l 



are all the terms that have resulted from dissociation. The corresponding pair 
of equations for [L n ] and [R n ] is automatically valid also for n = 1, so we 
have 

±[L n ] = ... + - (n - l) 7 s[L„], (17) 
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where 



^[Rn] = ... + - (n - \) ls [Rnl (18) 



N N N 

m=n+l m=n+2 m=l 



N N N 

W^=2 1S J2 [Rm]+lsJ2 [R m L 1 ]+ 1I [R n L 1 ]+SnlJ2i L mRl]- (20) 
m=n+l m=n+2 m=l 

Here, 5 n \ = 1 for n = 1, and <5 n i = for n > 2. 

We have calculated solutions using as initial condition [L n R] = 1 for 
different values of n and found that the evolution of [L n ] is very similar to 
that shown in Figure 1, so we do not need to reproduce this result here. 

3.3. Polymerization and dissociation 

We now add the polymerization equations of Sandars (2003) to Eqs. (13)- 
(18). Again, we begin by discussing first the homochiral case. In that case we 
have only two reactions, 

L n -i + L\ L n , (21) 
L n > L m -j- Z/ n _ m , (22) 

where is the reaction coefficient for attaching a monomer with the same 
handedness. The factor 2 on k$ indicates that polymerization can proceed on 
both ends of the polymer. This agrees with earlier approaches, and may be 
realistic for PNA polymerization, but not for RNA or DNA polymerization 
which usually proceeds only on one end. Since the monomer can be attached 
to any one of the two ends of the polymer, the overall reaction proceeds with 
the coefficient 2k s- The full reaction equations for the homochiral case can 
then be written as (for n > 3) 

= 2k s [L 1 ] ([L n _i] - [L n ]) + 2 1S ]T [L m ] - (n - l) 7s [L n ], (23) 

m=n+l 



while for n = 2 we have 

^1 = fc s [Lx] ([L X ] - 2[L 2 ]) + 2 1S E [L m ] - 75^2], (24) 

m=3 

where an extra 1/2 factor has occurred in front of the [L\] 2 term. (For n = 2, 
two pieces of the same species react with each other, whereas for n > 2 there 
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Figure 2. Isotactic equilibrium states with polymerization and dissociation, for different 
values of the universal parameter M changing over a range of six orders of magnitude. 

are always two different species, i.e. monomers and polymers; see BAHN for 
a more detailed discussion). For n = 1 we have 

, JV-l N 

-[In] = -2fc s [Li] i L n] + 275 Et L «]- ( 25 ) 

n=l n=2 

Note that this problem is governed by three parameters, ks, 7s, and the con- 
served quantity El, which only depends on the initial condition. These three 
parameters can be combined into a single non-dimensional parameter, 

M = ELks/^js (homochiral case), (26) 

that characterizes all possible solutions. Moreover, these equations possess a 
unique equilibrium state which is in general different for different values of 
A4; see Figure 2. Here we have normalized [L n \ in terms of js/ks to make 
it dimensionless. 

Given that there is a non-dimensional parameter (M) in the problem, there 
is no unique choice for a non-dimensional representation of time. Possible 
non-dimensional combinations are (as used in Figure 1) and Eikst. 
In Figure 3 we show the time dependence of [L4] (normalized by El) as 
a function of jst toward the equilibrium solution shown in Figure 2. Note 
that the approximate position of the maximum is always at around the same 
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Figure 3. Relaxation phase toward the isotactic equilibrium states in the presence of poly- 
merization and dissociation, for the same six different values of the parameter M as in 
Figure 2. 

value of 7st for values of M changing over six orders of magnitude. This 
shows that the typical relaxation time scale is governed by jg 1 . 

It is somewhat surprising that with dissociation, [L n ] always peaks at small 
values of n (at n = 2 or M > 10 or at n = 1 for smaller values of M). This 
can be quantified in terms of the mean polymer length Nl that can be defined 
as Nl = n [L n ] I J2 [L n ] (see BAHN). The resulting values of Nl approach 
3 for large values of M, but are otherwise always less than 3. 

3.4. Coupling to a substrate 

It is natural to proceed as in the model of Sandars (2003) and couple the 
polymerization equation to a substrate from which new monomers can be 
produced in a catalytic fashion. It is sufficient to discuss first the isotactic 



Table I. Mean polymer length Nl for different values 
of M, for isotactic polymers (here left-handed). 



M 


0.1 


1 


10 


10 2 


10 3 


10 4 


N L 


1.09 


1.55 


2.45 


2.92 


2.99 


3.00 



paper.tex; 9/02/2008; 3:22; p. 9 



10 



Brandenburg et al. 



case, so we add a source Ql to the right hand side of Eq. (25), 

JV-l N 

rr,j = Or - 9.ko\L-,] 

dt 



^[Li] =Ql~ 2k s [L 1 ] [Ln] + 2 7s £[L n ], (27) 

n=l n=2 



where Q l quantifies the source of new left-handed monomers. Since this term 
provides of source of left-handed building blocks, El is no longer conserved. 
Instead, as discussed by BAHN, El obeys the evolution equation 

^ = Q L -2k s [L 1 }[L N }. (28) 

In the absence of dissociation, a homochiral steady state is possible, where 
[L n \ is constant for all n > 2, so [Ln] is finite and is balanced by 

2k s [L 1 ][L N \. 

Obviously, Ql should depend on the concentration of the substrate, [S], 
so it is natural to write Ql = kc[S]CL, where Cl determines the efficiency 
of the production of left-handed monomers from the substrate. Since this 
generation is supposed to be a catalytic process, Cl should depend in some 
way on [L n ] itself; here we assume Cl = El, but different proposals have 
been made in the past (see BAHN for a discussion). The substrate itself obeys 
an evolution equation of the form 

^ = Q-(Ql + Qr), (29) 

where Q is a source for the substrate, and Qr in the present case. For the 
moment, this source can be thought of as being externally given, as in the 
model of Sandars (2003), but we will assume that this comes actually from 
the dissociation fragments. 

Regardless of the particular choice, we face a general problem in that dis- 
sociation causes the polymers to have finite length, so [Ln] ^ and hence no 
equilibrium state is possible any more. This causes a secular (linear) growth, 
so at some point the numerical integration develops an arithmetic overflow. 
An obvious way to balance this secular growth is to add a simple loss term, 
d[L n ]/dt = ... — 7[L„], where 7 is the degradation rate and the dots denote 
all the other terms that are already present. The result is show in Figure 4. 



3.5. Feeding the fragments back into the substrate 



Clearly, the dissociation model developed so far requires some modifications 
that are necessary to prevent the model from displaying secular growth when 
combined with the polymerization model of Sandars (2003) and to allow for 
an average polymers length of more than 3. One possibility would be to make 
the decay rate 75 dependent on n, for example in such a way that 75 = 
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Figure 4. Isotactic equilibrium states with polymerization, dissociation, and uniform degra- 
dation, for different values of M/W 4 (left), and the mean polymer length Nl (right), for 

7/7S = 20. 



for small values of n. One could also think of adding an overall loss term. 
Yet another possibility, that is close to our final proposal, is to recycle the 
monomers resulting from dissociation back into the substrate. In the end, 
however, we found it most plausible to assume that all fragments resulting 
from dissociation are recycled back into the achiral substrate. Thus, the source 
term would then be 

q = Wl + Wr + W lr + W rl + Wrlr + W LRL (30) 

where 

N N 

W L = Y, ^4 L) . W R = ]T nw™ , (31) 

n=l n=l 

is the total number of recycled building blocks (both left-handed and right- 
handed), 

N N 

W LR = (" + l Wn R) , W RL = J2(n+ lWn RL) (32) 

n=l n=l 

are the corresponding contributions from fragmented semi-spoiled polymers, 
and 

N N 

Wr L r = J2(n + 2) [R^ [L n R] , W LRL = £ (n + 2) [Li] [R n L] (33) 

n=2 n=2 

are the contributions from terminally spoiled chains. The new system of equa- 
tions is then 

±[L n ]=pW _( n _i) 75 [L n ], (34) 
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j t [Rn] = - (n - lhs[Rn], 05) 
^[L nJ R!] = p^) - {7/ + (n - 2) 75 } [L n i?x], (36) 

^[i? n Lx] = p^) - {7/ + (n - 2) 75 } [i^Lx], (37) 

where p^ , , p^ L ^ , and pl^ R ^ indicate the terms due to polymerization 
(see Sandars 2003, BAHN, Wattis & Coveney 2005); see the appendix. Equa- 
tions (34) and (35) are valid for all n > 1, but Eqs. (36) and (37) are only 
valid for n > 2. For n = 1 these equations reduce to 

^[L 1 R 1 ]=p[ LR) - 7 /[£ifli]. 08) 

We note that [i?iLi] = [L1JS1]. These equations are constructed in such a 
way that the total number (or mass) of right and left handed building blocks 
is conserved, i.e. 

M = [S] + E R + E L + E+ + E+ = const. (39) 

Here, E^ = J2( n ^ l)[^n] and E^ = J2( n ^ l)[-^n] have been introduced. 
We recall that due to recycling of the right and left handed building blocks 
through an achiral substrate, the total chirality, which involves E^ and E£, 
is not conserved; see Section 5 of BAHN. The quantity M can be expressed 
in non-dimensional form, 

M = Mkshs (general case), (40) 

which is conserved for all times. This is therefore the main control parameter 
of our model. It generalizes our earlier definition for the fully homochiral 
cases; see Eq. (26). 

In Figure 5 we show that increasing the value of M. leads to an increased 
range over which the racemic solution is unstable and a near-homochiral state 
emerges. Likewise, increasing the rate at which the spoiling monomers break 
off also broadens the range of permissible values of the fidelity for which the 
racemic solution is unstable. 

When / exceeds a critical value, the enantiomeric excess 

El~ E R 

V = (41) 

increases exponentially with time like e xt , where A is the growth rate. For 
/ = 1, the growth rate is (see Figure 6) A = 1.2A4 1 / 2 7s = (M/C575) 1 / 2 , so 
it is the geometrical mean between the polymerization rate Mks on the one 
hand and the dissociation rate 75 on the other. 
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Figure 5. The effects of M and 7/ on the bifurcation diagram. Increasing the values of M 
and 7/ allow near-homochiral states with decreased fidelity. 




Figure 6. Dependence of normalized growth rate X/'ys for the racemic solution and mean 
polymer length Nl for the homochiral solution on the total normalized mass parameter M. 

for f — 1, ki/ks = 1 and 77/75 = 1. 



Contrary to the model without recycling, the present model does allow for 
polymer lengths that can easily exceed the previous bound of 3. This requires 
large values of M; see Figure 6, where we plot the resulting values of Nl as 
a function of M. In fact, we find that to a good approximation, 

Nl or N R « 1.12M 1/4 (forM>10); (42) 

see Figure 6. We regard the possibility of long chains as a crucial property of 
any reasonable polymerization model. Furthermore, the fact that the model is 
now fully self-contained (M is conserved) makes it an appealing alternative 
to previous models. 
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4. Conclusions 

Dissociation of polymers appears to be an important component of any poly- 
merization model. The present work has shown, however, that the straight- 
forward usage of dissociation fragments for further polymerization does not 
yield realistic model behavior, because the maximum polymer length would 
not be more than 3. Various other modifications that could allow for longer 
polymers have been discussed, and it is likely that there are more possibilities. 
The main problem is that the fragments from dissociation tend to produce 
excessive amounts of short polymers that cause the average polymer length 
to be very short. Consequently, we have postulated that the fragments result- 
ing from polymerization are instead recycled into the substrate. The average 
polymer length then depends on the normalized dissociation time. In this 
model, no external source of the substrate is required, so the model is now 
fully self-contained. 

The model is governed by the total number of left and right handed ho- 
mochiral building blocks, the reaction rates for polymerization with the same 
and the opposite chirality, and the corresponding dissociation rates. These 
numbers can be combined into a single non-dimensional number that char- 
acterizes the behavior of the system. At the moment we have no clear idea 
about its value, but laboratory experiments should be able to determine not 
only this coefficient, but they should also allow us to test various aspects and 
predictions of the model. 

We recall that in order to draw conclusions about the time scale on which 
homochirality can be achieved, it is important to discuss the spatial extent of 
the system (Saito and Hyuga 2004b). Homochirality may develop rapidly at 
one point in space, but the handedness may be different at different locations. 
The relevant time scale for achieving global homochirality is therefore much 
longer and is given either by the diffusion time scale, which is very long, 
or by a turbulent turnover time which can be much shorter if turbulent flows 
are present (Brandenburg & Multamaki 2005). Obviously, the generalizations 
given in the present paper can directly be applied to their model provided the 
local growth remains still large enough. 

Appendix 

A. Polymerization terms 

In this appendix we state the terms describing the polymerization process. 
These terms are equivalent to those given and discussed in BAHN, Equa- 
tions (20)-(27). For n > 2 we have 

j4 l > = 2k s [L 1 ] (4 1 / 2 ) [L n _i] - [L n ]) - 2k I [L n ][R 1 ], (43) 
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pW = 2k s [R 1 ]( 


a (l/2) 


[Rn-l] ~ 


- [R n ]) - 2k I [R n ][L 1 ], 


(44) 


p n RL) 


= ks[Ri](a^[R n . 


-iL]~ 


[RnL]) 


+ ki[L 1 ](2[R n ]-[R n L]), 


(45) 


A LR) 


= ks[L 1 ](aW[L n . 


-iR] - 


[L n R]) 


+ k I [R 1 ](2[L n ] - [L n R]), 


(46) 


whereas for n = 1 we have 


(L) 
Pi = 


■ -Xl[L 


i] mdpi^ = —\r[Ri], where 



N-l N N-l N 

X L = 2k S ]T [L n ] + 2k 2 ]T [Rn] +ksY, i L nR] + ^/ E i R n L i ( 47 ) 
n=l n=l n=2 n=2 



N-l N N-l N 

X R = 2k s E [Rn] + 2kj E [Ln] + k s ]T [RnL] +k T J2 [LnR], (48) 

n=l n=l n=2 n=2 

and p\ RL ^ = p\ LR ^ = 0. 
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